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EXPLORING THE SPACE OF COMPACT SYMMETRIC CMC SURFACES 
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Abstract. We map out the moduli space of Lawson symmetric constant mean curvature surfaces 
in the 3-sphere of genus > I by flowing numerically from Delaunay tori with even lobe count via 
the generalized Whitham flow. 


1. Overview. In this note we map out a portion of the moduli space of embedded constant mean 
curvature (CMC) surfaces in the 3-sphere experimentally by a numerical implementation of the 
generalized Whitham flow [2] . This provides numerical evidence for the existence of the flow reaching 
arbitrary genus. 

Experiment. For each pair of integers g > 1 and n > 0 we construct numerically a 1-parameter 
family Eg of compact Alexandrov embedded CMC surfaces of genus g in with n controlling the 
lobe count: 

• the family starts at the Lawson surface and converges to a chain of two minimal 
spheres; 

• the family Eg (n > 1) converges at one end to a chain of {g + l)n CMC spheres and at the 
other to a chain of {g F l)n F 2 CMC spheres. 

Each surface in Eg has a cyclic symmetry of order g F 1 with four fixed points. 

These E'^ families were computed numerically via the generalized Whitham flow [2] , a topology¬ 
breaking flow through CMC surfaces in which starts at CMC tori and, as indicated by numerical 
evidence, reaches closed CMC surfaces of arbitrary genus. 

The generalized Whitham flow passes through each of the families Eg with n fixed and g in¬ 
creasing arbitrarily, starting at the tori Ef. To describe this initial data, recall [4] the embedded 
CMC tori in the 3-sphere consist of the 1-parameter family of homogeneous tori of increasing mean 
curvature starting at the minimal Clifford torus, along which bifurcate 1-parameter families of n- 
lobed Delaunay (equivariant) tori at sequencial bifurcation points Prn (see figure]^. The initial 
family E^ is made up of the homogenous tori between the Clifford torus and /32, together with the 
2-lobed Delaunay tori. The initial family Ef {n > 0) is made up of the (2n)-lobed Delaunay tori, 
the homogeneous tori between /32n and /32n+2, and the (2n + 2)-lobed Delaunay tori. 

The topology-breaking flow is described qualitatively as follows. The initial torus in Ef has a 
cyclic symmetry of order two with four fixed points. The flow preserves the topology of the torus 
minus two disks, formed by introducing two cuts connecting the fixed points in pairs. The flow 
retains the rotational symmetry, decreasing its angle a from tt to 0. When a = ^njiyg + 1), ^ G N, 
the surface can be completed by the rotational symmetry to a closed compact unbranched surface 
of genus g. At other angles a G 27rQ, the surface can be completed to a closed surface branched at 
four points. 

The generalized Whitham flow starting at a (2n)-lobed Delaunay torus can flow either to 
or E^. This flow direction is determined by the choice of order two symmetry: the symmetry with 
axes through necks flows to while the symmetry with axes through bulges flows to E'^ (see 

figure 1^. 
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Figure 1. The above plots show the mean curvature (left) and Willmore energy (right) along the 
family of genus 2 CMC surfaces depending on its conformal type encoded as the poleangle 
6 G [0,1/8]. For each plot the family starts at the Lawson surface (right) and ends at the 
chain of two spheres (left): the four poles of the DPW potential (umbilics of the surface) are 
^g±27ri0 events from right to left are marked by vertical lines: 1. Branchpoint on unit circle, 

2. Willmore energy = Stt, 3. Maximum mean curvature, 4. Maximum Willmore energy (> Stt). 



Figure 2. The moduli space of Lawson symmetric CMC surfaces in arising from even-lobed 
Delaunay tori. The vertical line represents homogeneous tori of increasing mean curvature starting 
at the minimal Clifford torus at the bottom. The horizontal lines represent Delaunay tori with 
increasing even number of lobes (bottom to top). The curved lines represent the 1-parameter 
families ^2 of genus 2 surfaces. The ffow occurs in the shaded regions; The families of genus 
g > 2 (not shown) are obtained by continuing the shaded regions beyond the curved lines. 
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2. The potential. We construct the families numerically via the generalized Weierstrass repre¬ 
sentation [T] for CMC surface in The Weierstrass data consists in a SI2C loop-valued potential 
^ with appropriate asymptotics in the spectral parameter A. The CMC immersion is obtained as 
F(Ao)T“^(Ai), where Aq, Ai G are the sym points^ and F is the unitary factor of the loop group 
Iwasawa factorization of $ solving the ODE d^ = 

The potential ^ for the families is a Fuchsian potential on CP^ with four simple poles 


(0.1) 


3 


c=E 

k=0 


Akdz 
z- Zk 


with order 2 symmetry (5*^ = 6{z) := —2), a := diag(i, —i) and real symmetry A) = 

The asymptotics of ^ in the spectral parameter A is determined by the requirements of the 
generalized Weierstrass representation and the Hopf differential of the surface: the upper right entry 
of the residue Aq has a simple pole at A = 0, the lower left entry of Aq has a simple zero at A = 0, 
and the residues of ^ have no other poles in the unit disk in Ca- 

To reach the families, we impose the condition that the eigenvalues Tuq, ±ui, be real and 
A-independent, with uq G (0, 1/4], ui G [1/4,/I), uq + = 1/2. This condition arises due to the 

fact that the eigenvalues control the angle a of the rotational symmetry being opened. With these 
assumptions, the potential ^ is a simpler replacement for the potential described in equation 2.1 
in [3] , to which it is gauge equivalent. 

We note that the points in the punctured unit disk in Ca at which the parabolic structure corre¬ 
sponding to ^ is unstable are those points for which the two eigenlines of Ai and ^42 corresponding 
to the positive eigenvalues coincide. We observed that the family ^ > 1 has n unstable points 
in the unit disk. 


2.1. Geometric parameters. The data determining a CMC surface in via the generalized Weier¬ 
strass representation is its potential two sym points in and the initial condition for the ODE 
d4> = 4>^. For Eg families, this data consists of three geometric parameters together with acces¬ 
sory parameters (coefficients of the residues of ^). The initial condition for the ODE d4> = is 
determined as the diagonal unitarizer of the monodromy, unique up to isometry of S^. 

The three real geometric parameters (7, c, H) are 

• the angle a := Attuq of the rotational symmetry being opened; 

• the conformal type 7 := [zq, — Zi, —Zq, Zi] G M of the four punctured CP^, 

• the mean curvature H := i(Ao + Ai)/(Ao — Ai). 

2.2. Accessory parameters. The eigenvalues of the residues of determining the angle of the rota¬ 
tional symmetry, must be controlled during the flow. Since the complex dimension of the space of 
monodromy representations on the 4-punctured sphere with fixed eigenvalues is, roughly speaking, 
2, the residues of ^ can be parametrized by two meromorphic functions x and ^ of A. For numerical 
calculations, we represent x and y as truncated power series in A at A = 0. Because the monodromy 
of ^ is to be evaluated on the unit circle we require that the potential ^ is holomorphic in 
the punctured unit disk. This holomorphicity is achieved by the introduction of polynomials and 
constraints on these polynomials. 

More precisely, let x, ^ be functions on the unit disk with f, 1 /t and y holomorphic. Let p/c, 

(/c G {0,..., 3}) be polynomials satisfying the constraints that pk monic and 

(0.2) Gjk ■■= {PjQj - Vj) - {PkQk - Vk) , (i, fc e {0, ...,3}) 

vanishes. Under this constraint, the residues of the potential ^ can be parametrized by p/^, g/c,T, p, 
(A: G {0,..., 3}) as 

y -2/12/3 A 
PiPsx y 

y := Vk+PkVk , 

where U 2 := —uq and U3 := —ui. 


(0.3) 

(0.4) 


Aq = 


Ai — 


-y P0P2/A 

.-2/02/2A y 
P ■= P0P1P2P3 , 2/fe := Qk + py/Pk , 
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Figure 3. Genus 2 surface in with five lobes (left). Genus 2 surface in with six lobes (right). 


For numerical computation, the accessory parameters are 
(0.5) A := (coeffpo ,... coeSps | coeff go, • • •, coefF gs j xq, ■ ■ ■, xiv I yo, ■ ■ ■, Vn) 

where coefF q denotes the coefficients oF a polynomial g, and the series 

OO OO 

(0.6) ^ = X] 

k=0 k=0 

are truncated to power N. The constraints ( |Q.2[ ) are 

(0.7) Ca •= (coeff eoi, coeff eo 2 , coeff eos) . 


3. The flow. The generalized Whit ham flow is defined to preserve intrinsic and extrinsic closing 
conditions. This flow is an implicit infinite dimensional ODE computed numerically by truncation 
to a finite implicit system of the form AX + 5 = 0; X is obtained as the least squares solution to 
this system. The coefficients of the system depend on the monodromy of the potential computed 
by a separate nested ODE. 


3.1. Intrinsic closing condition. The intrinsic closing condition is that the monodromy of ^ is uni- 
tarizable along the unit circle More concretely, let {k G {0, ...,3}) be the monodromy 
generators for ^ along a curve based at 2 ) = 0 which winds once counterclockwise around the pole Zk, 
and let tjk •= | tr By proposition 2 in [3], the monodromy is unitarizable when tij G (—1, 1). 

Hence we impose the constraint along 


( 0 . 8 ) 


cj := (Imtoi ,Imto 2 ,Imto3 ,Imti3 ,Imt23) • 


Eor numerical computation, this constraint is implemented by imposing the condition (0.8) at S 
equidistant sample points ii^ — {k G {0,... S' — 1}) along S^. Eor the flow, the number S of 

sample points must be large relative to the number N. This is the vanishing of 


(0.9) 


Cl := (c/(/io), • • •, c/(/is_i)) . 
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3.2. Extrinsic closing condition. The extrinsic closing conditions are that every monodromy of the 
unitary frame M satisfies M(Ao) = M(Ai) G {± 1} at the two sympoints Aq, Ai G S^. By proposi¬ 
tion 1 in [3] , this is the condition that 

(0.10) ce := [4, h, ^2, 4] - [^0, -^1, -^0, ^i] 

vanishes to first order in A at each of the two sympoints, where ik ^ {k G {0,...,3}) is the 
eigenline of Ai^ corresponding to its positive eigenvalue. With prime denoting the derivative with 
respect to A, this is the vanishing of 

(0.11) Ce := (ce(Ao), c'^{Xo), ce{Xi), . 


3.3. The flow. The flow is defined in terms of the function /, defined to vanish when the intrinsic 
and extrinsic closing conditions are satisfied: 

intrinsic closing condition Cj 
extrinsic closing condition Ce 
accessory parameter constraint Ca 

To flow along a curve in the 2-dimensional isosurface /“^(O) we consider 
(0.13) t A (t, u, ( 7 , a, H, A) A {Ci, Ce, Ca) , 

where t is the real flow parameter, r is a real free parameter, h is an explicit immersion controlling 
the direction and speed of the flow, and Y is the sought function defined implicitly by the condition 
F oY = 0, where F := f o h. With dot denoting the derivative with respect to t, Y is defined by 
the implicit ODE 

(0.14) dFY = 0. 


( 0 . 12 ) 


geometric parameters (7, c, H) 
accessory parameters A 


In matrix form, 
(0.15) 


[B A] 


= 0 , 


that is AX + 5 = 0. 


The vector field X is obtained as the least squares solution to the system AX + 5 = 0. 

The generalized Whitham flow starts with the initial data for a homogeneous or Delaunay torus [2] , 
with (7, c, H) = (constant, t, r), reaching a surface in by increasing genus. Starting from such 
a surface, the Whitham flow moves along the Eg family, with (7, c, H) = (t, constant, u). 

4. Lawson symmetric surfaces. A Lawson symmetric CMC surface is a compact CMC surface in 
of genus g > 1 which enjoys a cyclic symmetry of order g 1 with four fixed points. The 
families Eg described above are Lawson symmetric, and have an additional symmetry induced by 
the hyperelliptic involution. We conjecture, for each pair of integers g > I and n > 1, the existence of 
an additional 1-parameter family Eg of Alexandrov embedded Lawson symmetric surfaces which lack 
the symmetry induced by the hyperelliptic involution. This family is reachable via the generalized 
Whitham flow from the (2n + l)-lobed Delaunay tori, and converges to a chain of {g + l)n + 1 CMC 
spheres. 

Conjecture. The space of Alexandrov embedded Lawson symmetric CMC surfaces consists of the 
families Eg and Eg. 

In the case g = I, this moduli space is connected. In the case of fixed g > the families Eg 
ranging over n G N are disconnected from each other. 
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